A 1/t damped electrostatic electron plasma wave 
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In an electron plasma, besides the exponentially Landau damped electron plasma waves, there 
can also exist modes, namely the ballistic or Case- Van Kampen modes, that usually decay faster 
(say, ~ exp(— t 2 )) than the linear Landau damping. In this paper, a slower (namely ~ 1/t) damped 
mode is considered. The latter resolves the paradox of why Landau damping is difficult to realize 
in the numerical simulations based solely on the Vlasov-Ampere (V-A) equations, even though the 
modes have the same dispersion properties as that obtained from the Vlasov-Poisson equations. The 
mode of interest here corresponds to a residual mode of the V-A system, and its spectrum is also 
discussed. 

PACS numbers: 52.35.Fp, 52.25.Dg, 52.65.-y, 52.30.-q 



I. INTRODUCTION 

The evolution of linear electron plasma wave pertur- 
bations in an unmagnetized plasma is generally deter- 
mined by an initial-condition independent part (the nor- 
mal mode) governed by a dispersion relation and an 
initial-condition dependent, or transient, part. The for- 
mer is Landau damped [T] and decays exponentially. Van 
Kampen[5] and Case[3] showed that there can exist other 
modes, namely the ballistic or Case- Van Kampen modes. 
These are undamped in the sense that there are no solu- 
tion of the linear Vlasov-Poisson equation with 9(w) < 
for the harmonic perturbation ~ cxp(— iwt), where w is 
the mode frequency. However, because of phase mixing 
(or the Riemann-Lebesgue lemma) , the perturbations are 
still damped when all such modes are considered. The 
Case- Van Kampen modes have been investigated numer- 
ically and analytically, and interesting and important re- 
sults have been foundjIHH]. For example, Ng et o/.[H[5] 
found that the continuous Case- Van Kampen spectrum 
can disappear and be replaced by a discrete spectrum 
close to the Landau solution. Recently, Bratanov et 
al. [51 [7] reconsidered the problem and new insights on the 
relation between the Landau and the Case- Van Kampen 
modes were found. 

The ballistic modes usually decay faster than Lan- 
dau damping, especially when the initial perturbation 
is an entire function [^ 110] . For example, for the ini- 
tial perturbation Sfa = AoCxp[— [u — u a ) 2 /u1], where 
\u a \ ^> (l u fcUw r /fc|), the ballistic mode decays like 



exp(— iku a t 



It 2 ) [TTJ. As a result, the asymptotic 



behavior of the plasma is determined by the Landau 
damped normal modes. The ballistic modes are thus ig- 
nored in most studies on plasma waves. Since they are 
short lived, the ballistic modes are also difficult to find 
in the experiments. However, there can exist slowly (rel- 
ative to that of Landau) damped ballistic modes if the 
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initial perturbations are not entire functions. 

In this paper, we are interested in an algebraic decay 
mode excited by the initial perturbation 



Sfo = 



A , for \u- u a \ < u b , 
0, otherwise, 



(1) 



which is not an entire function. Asymptotically the cor- 
responding ballistic mode ([I]) behaves like [TT] 



oc - sm(kubt) exp(— iku a t), 



(2) 



so that it decays like 1/t. Our main purpose here is follow 
the evolution of such a slowly decaying ballistic mode. 
Since it is difficult to obtain the asymptotic behavior with 
the standard eigenvalue formulation, we shall consider 
the initial value problem. 

In numerical simulations of electrostatic waves in plas- 
mas, the Vlasov-Ampere (V-A) equations are often more 
convenient to use than the Vlasov-Poisson (V-P) equa- 
tions since the Poisson's equation is elliptic. However, 
when considering electrostatic one-dimensional (ES1D) 
problems in Vlasov plasmas, the results from the V-A 
system may violate the Gauss's law. To overcome this 
problem in their numerical simulation of Landau damp- 
ing, Home and Freeman [T^] used the V-P equations for 
the first few steps of the numerical integration, before 
continuing with the V-A equations. Here we show that 
the V-A system alone, i.e., without the Home- Freeman 
Poisson start, can also yield solutions that are equivalent 
to that from the V-P system. Our result also explains 
why it is difficult to obtain Landau damping using the 
V-A system. 



(a)k=0.5, dv=0.15625 



(b) t=99.99, dt=0.01 




Reshape 8f(k,v) to f(x,v)=f (v)+5f(x,v) 
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FIG. 1: Benchmark of the Landau damping solution. 
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FIG. 2: Simulation of the 1/t decay ballistic mode. 



II. SIMULATION OF THE 1/i-DECAY 
BALLISTIC MODE 



The normalized linear Vlasov-Poisson system is given 



dt&f = -ikvSf + SEd v f , 
ikSE = - / Sfdv. 



(3) 



where /o = exp(— v 2 /2)/^/2tt. 

We solve pi) as an initial value problem from t = 
to tend = NtAt using a 4th-order Runge-Kutta scheme. 
The discrete velocity space is from t> m ; n to w max - There 
are N v uniform grids of size Av — (u max — v m - ln )/N v . 



> 




FIG. 3: Reshape 8f(k,v) to f(x,v) of the Fig. [5] ballistic 
simulation result. 



As mentioned, for most initial perturbations Sf(t = 0) 
or 5E(t = 0), the asymptotic behavior of 5E(t) is de- 
termined by the normal mode that is Landau damped. 
We shall first verify this solution. For, k = 0.5, the 
least damped Landau solution is Wtheory = co r + ij = 
1.4157 — 0.1534*. We use as initial perturbation a Gaus- 
sian Sf with A = 0.05, u a — 1.0, and Ub = 1.0. The 
electric field 8E(t = 0) and 8E(t) can be obtained from 
the Poisson equation using the initially given and the 
calculated Sf, respectively. The simulation results are 
shown in Fig. [I] We can see in the panel (d) that 
the simulated curve for Landau damping matches well 
with the theory. With dt — 0.01 accuracy, we have 
w s i mu iation = 1.42 — 0.15*. A non-physical recurrence ef- 
fect (the Poincare recurrence) can also be seen in the 
panels (c) and (d) at Tr = 2n/kAv, which is due to the 
discreteness of the velocity space [T3] ■ 

The 1/t damped ballistic mode from the initial pertur- 
bation (fTl) can be seen in Fig. [2] for ^o = 0.05, u a = 0.0, 
and Ub — 2.0. As reference, the red dashed line in the 
panel (d) is for y — 3.3a;. We see that \5E\ decays as c/t, 
as analytically predicted in (pi). Note that the constant 
c depends on Aq. To avoid the non-physical recurrence 
effect, we have used a very small Av in this case. 

In the above we have obtained the 1/t damped ballistic 
mode from the V-P system. Our linear simulation was 
carried out in the fc-space. For completeness, we give 
in Fig. p^the corresponding (x,v) phase space by direct 
mapping of Fig. fusing the relation f(x,v) = fo(v) + 
$l[6f(k,vy kx }. 



III. ANSWER TO THE ES1D 
VLASOV-AMPERE PARADOX 

We now consider the ES1D Vlasov-Ampere paradox 
described in the Introduction. 



A. E modes 

The Poisson's equation in ([3| is now replaced by the 
Ampere's law 



(a) k=0.8, dv=0.0097656 



(b) t=24.99, dt=0.01 



d t SE = / vSfdv 



(4) 



From the charge continuity equation dtp + d x J = 0, it 
is easy to show that the V-A and V-P systems are equiv- 
alent if the Poisson's equation is used to supplement the 
initial condition in the V-A system, as discussed above. 
The dispersion relation from the V-A system is 



Dy A (w,fe) 



vdfg/dv 

uj — kv 



dv = 0, 



and that from the for the V-P system is 

dfo/dv 



Dyp(u,k) = k 



kv 



dv = 0. 



(5) 



(6) 



It is readily shown that ((5| and ^ are equivalent, and 
they should thus both yield the normal modes given by 
the dispersion relation, which is independent of the ini- 
tial condition. That is, one should also be able to get 
the Landau damped solutions from the V-A equations 
without using the Poisson start. However, existing sim- 
ulations indicate that it is difficult [T2] . 

In view of the above discussion, we now show that the 
ballistic mode (here we mean modes related to the ini- 
tial perturbation, not necessarily that of the traditional 
definition ui = ku) is the key factor in this paradox, i.e., 
the asymptotic behavior of the V-A system is mainly de- 
termined by the ballistic mode, instead of by the normal 
mode as for the V-P system. 

Applying Laplace transform (L p ) in time and Fourier 
transform (F r ) in space to the V-A system, we can get 



SE(t,k) = L- l 6E(oj,k) 



(7) 



t doj 



c u 



2tt 



du 



uSf k (0) 
(w — ku)D(uo, k) 



and for the V-P system we get 

f e- iult du 



6E(t, k) = 



2vr 



du 



Sf k (P) 



(w — ku)D(uj, k) 



SE k (0) 
D{uj,k) 



(8) 




simtl&tb? tpfte[8 E] 
simulatbn lh|lm[8E][ 
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FIG. 4: V-A simulation of the residual E mode 
(a) k=0.8, dv=0.0097656 (b) t=24.99, dt=0.01 




FIG. 5: V-A simulation of Landau damping using Poisson 
start 



We are interested in the asymptotic behavior, so that 
only the to = ku pole and the maximum-Sw m normal 
mode are kept. From fff\ we have, 



From ([8]) we have 
lim SE = — 

t—$-oo 27T 



(111: 



„ — ikut 



Sfk(0) { er^mt 

Dyp(ku, k) {u m - ku)d u 



D 



ballistic part 



lim SE = 



1 

2^ 



du 



ikul. 



Sf k (0) 



DypJJzu, k) (w m — ku)d Um D 



Landau part 

(10) 

wDyp , we can see that the 

main difference between (111) and (jlOh is in the initial 



Noting the relation kDyj^ 



ballistic part 
e~ wt duj SE k (0) 

v ' 

initial 5E part 



Landau part 



5E, or ^(O), which is a constant. To see the typical 
asymptotic behavior, we let Dyj^ s=a u in the 5E part of 
the integral and obtain the residual mode 

(9) Hm [^ VA (t)-^ VP (t)+5£p oisson (0)] ex SE(0). (11) 
On the other hand, if we treat -DyA as constant, we 
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FIG. 6: V-A simulation of the E mode using initial perturba- 
tion (fl2b 



(a) V-P eigenmode solutions, k=0.8 < b > V ~ p s P ectral densi, y' v max= 8 
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(c) V-A eigenmode solutions 
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FIG. 7: Eigenmodes in V-A and the V-P systems with the 
comparison to the Landau damping normal modes 



will get lim t _ > . 00 ((5£ , Y^ — SE\rp) oc 1/t. In both cases, 
Landau damping is not evident. Similar to the ballistic 
mode, this mode is also from the initial perturbation and 
can be considered as a Case- Van Kampen mode. We 
shall call it the E mode, in order to distinguish it from 
ballistic mode. 



Using ( 11 ), one can separate the Landau and the resid- 



ual E modes via a novel setup and diagnostics strategy. 
An example of separating the Landau and E modes is 
given in Fig. [I] for SE(0) = 0.02 and the same initial 
Sf perturbation as in Fig. [I] In the panels (c) and (d), 
the real and imaginary parts of SE can be seen as the 
Landau and E modes, respectively. If SE(0) is purely 
imaginary, the opposite holds. However, in practical sit- 
uations it may be difficult to separate the Landau and 
E modes, especially in experiments or nonlinear simula- 
tions. Moreover, other profiles of (5/(0) or SE(0) can also 
lead to different results. For example, if u a = instead 
of u a = 1.0, the Landau damping in the above simulation 
would vanish. That is, as expected, the result is sensitive 
to the initial perturbations. 

The results for a V-A simulation of Landau damping 
using a Poisson start is given in Fig. [5j The parameters 
are the same as that for Fig. |4| except that a V-P run is 
added in the initial i t time step(s). In this linear simu- 
lation, i t = 1 is sufficient. For nonlinear simulations one 
may need more time steps [12]. In Fig. [5l we can see 
that both the real and imaginary parts oiSE match the 
Landau result. 

In the V-A simulation, we found that the residual E 
mode is dominant. The initial perturbation (nj) also can- 
not prevent it. However, we found that the perturbation 



can result in a roughly-1/i decay mode when u a = 0, as 
shown in Fig. [6] for u\, = 5.0. 



B. Eigenmodes in the V-A system 

There are also other phenomena related to Landau 
damping. For example, it is known that the Landau 
damped normal mode is not an eigenmode in the V-P sys- 
tem, although a growing normal mode can be an eigen- 
mode (see, e.g., Ref. 3J. BratanovJSj recently reported 
new numerical investigations of the connection between 
the Case- Van Kampen eigenmode and the Landau nor- 
mal mode. It is found that there is a spectral density 
accumulation around the real frequency of the Landau- 
damped mode. 

If we rewrite the governing equations into the ma- 
trix form M ■ F = luF, with dt = —iuj. The eigen- 
vector F for the V-P system is F = {Fj} = {Sf(vj)}, 
where Vj = v m i n + (j — 1)A« and j — 1, 2, 3, ...,N V + 1. 
Bratanov[6] found that although all the eigenvalue solu- 
tions of F are undamped, the integral of them, which 
is related to SE, can yield Landau damping because of 
phase mixing, similar to the Van Kampen modes. This 
can also be seen in the panels (a) and (b) in Fig. [7] 

The question is then: what if SE is also contained in 
the eigenvector Fl Will it give the Landau solution di- 
rectly, so that the Landau damped mode can also be an 
eigenmode? The V-A system can combine SE into F di- 
rectly. The eigenvector of the V-A system then becomes 



,N + l, 



Sf(0) = 



A u, for \u-u a \ < u b , , 12 s 

0, otherwise, ^ ' 



p_ rp, _ f SfiVj), j = 1,2,3," -,., , ., 

* ~ ! jl ~\5E, j = N + 2. [L6) 

Using pi and Q, we can easily find the elements of 



the eigenvalue matrices M for the V-A and V-P systems, 
respectively. The solutions are shown in Fig. [71 However, 
we see that there is no eigenmode, and the normal modes 
are overlapping. That is, our result also indicates that 
the Landau damped mode in the V-A system is not an 
eigenmode. Apparently, the only difference between the 
V-A and V-P systems is in the spectral density around 
uj r = 0: there is an extra accumulating point at u r = of 
the V-A system [compare the panels (b) and (d)]. Noting 
(111, we can attribute this new accumulating point to the 



E mode, in particular the constant residual E mode. 



and nonlinear effects are precluded. Further investiga- 
tion (say, using particle-in-cell simulation or kinetic the- 
ory) including these effects as well as a more realistic 
initial pulse, are still needed in order to see if such a 
mode can indeed be produced experimentally. We have 
also resolved the ES1D V-A paradox and found an inter- 
esting mode, namely the residual E mode, which does not 
seem to have been investigated earlier. However, we note 
that, when considering the vibrations of a plasma in an 
external eletric field, Landaup] found analogous results 
as that in the V-A system. 



IV. SUMMARY 
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